Random Matrix Theory and Classical Statistical Mechanics: 

Spin Models 

H. MeyerB and J.-C. Angles d'Auriacfl 

Centre de Recherches sur les Tres Basses Temperatures, 
BP 166, 38042 Grenoble, France 
(PREPRINT December 18, 1996) 

We present a statistical analysis of spectra of transfer matrices of classi- 
cal lattice spin models; this continues the work on the eight-vertex model of 
the preceding paper. We show that the statistical properties of these spec- 
tra can serve as a criterion of integrability. It provides also an operational 
numerical method to locate integrable varieties. In particular, we distin- 
guish the notions of integrability and criticality considering the two examples 
of the three-dimensional Ising critical point and the two-dimensional three- 
state Potts critical point. For complex spectra which appear frequently in 
the context of transfer matrices, we show that the notion of independence of 
eigenvalues for integrable models still holds. 
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I. INTRODUCTION 

Random Matrix Theory (RMT) has been applied in surprisingly many fields of physics 
and mathematics. We have proposed in a recent letter [Q] its application to transfer matrices 
of lattice models in classical statistical mechanics. In a preceding paper ||, referred to 
hereafter as Paper I, we have given the details of this RMT analysis applied to the general 
eight-vertex model. We have numerically shown that the integrability of the model can be 
seen on the statistical properties of the entire spectrum of transfer matrices of vertex models. 
Using this as a criterion for integrability, we have located all the known integrable varieties 
in the parameter space. In this paper we continue this work with the study of spin models. 
Many aspects have already been presented in paper I, so we will recall below only the basic 
ideas of the RMT analysis with emphasis on the points which are specific to spin models. 

An important area of application of RMT is the characterization of chaos @,f|]. One 
can describe the fluctuations of energy spectra of chaotic systems with some ensembles of 
RMT, while the spectra of regular systems show the characteristics of independent numbers 
(Poissonian ensemble). In classical (Hamiltonian) mechanics, the notions of regular and 
chaotic dynamics coincide with the notions of integrability and non-integrability. But in 
quantum mechanics the notions of chaos and integrability are less precise; one nowadays 
adopts the criterion of RMT as a definition of quantum chaos. For models of quantum 
statistical mechanics one can adopt a definition of integrability related to the Bethe ansatz: 
an integrable system is a system for which a complete set of eigenstates having the Bethe 
ansatz form exists. For a classical statistical mechanics model, the notion of integrability 
is generally related to the Yang-Baxter equations. For example, solving the Yang-Baxter 
equations for the symmetric eight-vertex model (also called the Baxter model) allows to 
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build a one-parameter family of commuting transfer matrices, and finally to compute the 
free energy of the model ||||. The spin version of the Yang-Baxter equations are the 
star-triangle equations which have been introduced earlier Solving these star-triangle 

equations also allows to construct an infinite family of commuting transfer matrices. We 
note that a vertex model can be turned into a spin model with many-spin interactions (IRF 
model) (see for example [§]). In this paper, we treat spin models including only two-spin 
interactions, possibly coupled to an external field, and by integrable we mean Yang-Baxter 
integrable as well as star-triangle integrable. 

The spectrum of an integrable system, after a suitable treatment, has been shown to have 
many properties of a set of random independent numbers (Poissonian behavior), while the 
spectrum of a chaotic system is described quite accurately by the spectrum of matrices of 
statistical ensembles. The choice of the proper ensemble depends on very general symmetry 
properties of the model under consideration. For a time reversal symmetric model this 
ensemble is the Gaussian Orthogonal Ensemble (GOE) This classification scheme 

has been applied succesfully to quantum spin models in one dimension [|10| , |lTf and on two- 
dimensional lattices fl2|-|r5|j. 

For a classical lattice spin model, the energy spectrum is usually very simple. (For an 
Ising model this spectrum contains all the possible numbers of violated bonds i.e. the set 
of integers [0, iVj where Nf, is the number of bonds; the physical properties of the model 
are contained in the degeneracies, its statistical properties have been studied in Ref. [pf] .) 
Therefore the analysis has to be performed on another quantity. The transfer matrix is such 
an operator (related to the Hamiltonian) which describes completely the thermodynamic 



properties of a system including the size effects ||16||, and we will perform the RMT analysis 



on its spectrum. It has already been shown in paper I and in [|T| that the notion of Yang- 
Baxter integrability coincides with a Poissonian spectrum, and that nonintegrable spectra 
are described by GOE matrix spectra. To perform this RMT analysis in its usual form, 
one needs to deal with a real spectrum. For spin models and when the interactions are 
nonchiral (symmetric), this can be achieved using the so-called row-to-row (or layer-to- layer 
in higher dimension) transfer matrix as explained below. We also have studied the spectra in 
the nonphysical regime where the Boltzmann weights are not positive. Indeed, most of the 
analytical results concern varieties in the entire parameter space including the region where 
the Boltzmann weights are negative. In this region, even row-to-row transfer matrices can 
have complex eigenvalues. The necessary changes in the analysis will be briefly discussed. 

The plan of this paper is the following. In Sec. |TJ we recall the numerical methods of the 
RMT analysis with a special emphasis on the specificity of the spin models we investigate. 
In Sec. |n| we present the numerical results, succesively of the two-dimensional Ising model 
in absence of magnetic field (a paradigm of an integrable system!), the same two-dimensional 
Ising model in a magnetic field, the three-dimensional ising model (a paradigm of noninte- 
grable system!), and the three-state Potts model. The three-state Potts model provides an 
example of a point which is integrable and critical at the same time. By contrast, the three- 
dimensional Ising model provides a nonintegrable, but critical, point. Finally, we investigate 
a nonphysical self-dual point of the three-state Potts model. This is a first attempt to study 
complex spectra in the context of this statistical analysis of transfer matrices. We conclude 
in Sec. |V| with a discussion. 
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II. NUMERICAL METHODS OF RMT IN THE CONTEXT OF SPIN MODELS 



The machinery of Random Matrix Theory has been explained in details in Paper I. It 
consists in five distinct steps: 

• the first one is to choose a representation basis for the operator and to construct the 
corresponding matrix; 

• the second step consists in finding the parameter-independent stable subspaces and 
find the matrices representing the operator in each of these subspaces; 

• the third step is to diagonalize each matrix; 

• the fourth step is to 'unfold' each spectrum; 

• and in the fifth step all the spectral quantities are computed. 

In this section, we briefly discuss these five points and give the details which are specific to 
the spin models studied in this paper. 

(i) We use the transfer matrix formalism, where the lattice is built up adding identical 
"generating sub-lattices". These generating sublattices can be rows for two-dimensional 
models or rectangular layers for three-dimensional lattices. For spin models, it is well known 
that this transfer matrix can be factorized as: 

T(K U K 2 ,..., K n ) = V 1/2 (K 2 , . . . , K n )H(K 1 )V 1 / 2 (K 2 , ...,K n ), (1) 

where H contains the interactions between two generating sub-lattices and V contains the 
interactions inside a generating sub-lattice; V is usually a diagonal matrix. Ki are the 
coupling constants in the different directions, we assign the direction 1 to be the direction 
in which the lattice grows. For each case we study in this paper, the precise form of V and 
H are given in the corresponding section. In the form (0) the transfer matrix is a symmetric 
matrix if the interaction between the two layers and thus the matrix H is symmetric. If 
the Boltzmann weights are real positive, the entries of the transfer matrix are also real 
positive. In this case of a real and symmetric (transfer) matrix the spectrum is real so that 
the methods for the statistical analysis presented in paper I apply. 

(ii) To find the parameter independent subspaces one needs to know the symmetries of 
the generating sub-lattices. For the transfer matrix of a <i-dimensional hyper-cubic lattice 
these generating sub-lattices are (d — l)-dimensional lattices. For the square lattice, the 
generating sub-lattice is a periodic chain. As explained in Paper I the symmetry group is 
the dihedral group V^. For a three-dimensional cubic lattice we use a square lattice as 
the generating sub-lattice. For an isotropic square lattice the automorphy group has been 
detailed in |L5] (with emphasis on the atypical case N = 4). It is a large group which 
leads to consequent size reduction (see Sec. [111 C|) . In addition to these space symmetries 
there is a 'color symmetry' when no field is applied. This symmetry reflects the fact that 
the Hamiltonian is invariant under a permutation of the possible values of the spin variables 
(i.e. spin reversal symmetry for the Ising model). For a g-state model, this is a S q symmetry; 
in this paper we work with q = 2 and q = 3. This 'color-symmetry' commutes with the 
space symmetries. 
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(iii) The diagonalization of the blocks is done numerically using standard procedures 
of the LAPACK library. We note that the block-diagonalization of step (ii) requires more 
numerical effort (CPU time) than the diagonalization of one block. 

(iv) The unfolding procedure is the same as the procedure detailed in paper I. It produces 
the unfolded eigenvalues e from the raw eigenvalues A. The spectra of spin and vertex 
models are quite similar. However, we also have analyzed complex spectra. In that case the 
eigenvalues are seen as points in the plane (and not on a line) the local density of which 
has to be made constant. To unfold complex spectra we follow the procedure described in 
Chap. 8.6 of Ref. |§. 

(v) The spectral analysis is performed on the same quantities as in paper I. These are the 
level spacing distribution P(s) of the differences between two consecutive unfolded eigen- 
values Si = Cj+i — ti (for a nonintegrable model, the eigenvalue spacing distribution is very 
close to the Wigner surmise for the GOE: 

P(s = X i+1 - Xi) = exp(-7rs 2 /4) , (2) 

in contrast, with the exponential P(s) = e~ s for a set of independent eigenvalues for an 
integrable model); the spectral rigidity: 

Ag(L) = (-min/ (N u (e) - ae - bf ale ) , (3) 

where N u (e) = J2i^( e — £ i) is the integrated density of unfolded eigenvalues and (...)« 
denotes an average over a; another quantity of interest is the number variance S 2 (L) defined 
as the variance of the number of unfolded eigenvalues in an interval of length L: 

,2\ 



S 2 (L) 



N u (e + \) - N u (e - £ 



(4) 



where the brackets denote an averaging over e. The expected behavior in the limiting cases 
of independent numbers and of GOE spectra has been recalled in paper I; for more details 
see for example Refs. 

We also recall the parametrized probability distribution we use to interpolate between 
the Poisson law and the Wigner law: ch. 16.8]: 

p p ( s ) = c (l + 0) s p exp (-cs p+1 ) (5) 

with c 



. The interpolation parameter f3 proved itself to be a useful indicator 
for the localization of integrable varieties 0. There should be no confusion of this parameter 
(3 with the inverse temperature 1/ksT. 

III. RESULTS OF THE RMT ANALYSIS 
A. Two-Dimensional Ising Model without Magnetic Field 

We start our analysis with the two-dimensional Ising model, which is well known to be 



integrable in the absence of a magnetic field [ 18 1 . To work with symmetric matrices we use 
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the row-to-row transfer matrix. For a rectangular N x M lattice with periodic boundary 
conditions the (reduced) Hamiltonian reads: 



N-l M-l 

{3H° NM (K U # 2 ) = - E E °ii (KiVij+i + K 2 a i+ltj ) , (6) 

i=0 j=0 

where f3 is the inverse temperature, and K\ and K 2 are the coupling constants in the two 
directions divided by the temperature. The partition function is 

Z NM (Ki,K 2 ) = ]Texp (pH%, M {Ki,K 2 ) 

= TtT^(K u K 2 ), (7) 

where the row-to-row transfer matrix T^(Ki, K 2 ) is given by: 

T N (K 1 ,K 2 ) = V(K 2 )H(K 1 ) , (8) 

The matrix V(K 2 ) is a 2 N x 2 N diagonal matrix with entries 

N-l 

[v(K 2 )u >a = n ™r* +i , (9) 

i=0 

where aj = ±1 according to the i th digit of the binary representation of a and w 2 = e Kz . 
The matrix H(Ki) is a symmetric matrix 

h( Ki ) = h( Kl r N , h( Kl ) = ^ j_\ ^ j , (io) 

with Wi = e Kl . Note that K\ and do not play the same role here. If all the Boltzmann 
weights are real and positive, and using the circular property of the trace, T N (Ki, K 2 ) can 
be replaced by the similar matrix 

T N (K U K 2 ) = V 1 l 2 (K 2 )H(Ki)V 1 l 2 (K 2 ) . (11) 

It is clear that Tn(Ki, K 2 ) is symmetric and therefore has a real spectrum. This spectrum has 
been completely worked out for even size N in From the explicit form of the eigenvalues 
it is easy to check that the entire spectrum is invariant by negating the Boltzmann weights 
Wi and we can restrict ourselves to the physical case Wi > 0. 

In Ref. |19j it is shown that the problem of diagonalizing T^{K\, K 2 ) can be turned into 
a problem of free fermions. Namely one has: 

T N {K 1} K 2 ) = (2sinh2K 1 ) Af/2 exp (- £e g - ~) j , (12) 

where £ q are fermionic operators. The dispersion relation e q (Ki,K 2 ) is a cosine and its 
detailed form depends on the parity of the number of ^-particles. Therefore, fixing all 
the "quantum number" q i: where is the momentum of the i th ^-particle, completely de- 
termines an eigenvalue and an eigenstate of the transfer matrix. However, instead of using 
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the complete set of quantum numbers, we take only into account the quantum numbers 
which correspond to parameter-independent symmetries. As mentioned in Sec. 0, the space 
symmetries of T N (Ki, K 2 ) form a group isomorphic to V N = Z N txZ 2 , and in the absence of 
a magnetic field the spin reversal symmetry induces an extra Z 2 symmetry which commutes 
with T>n. Using the projectors independent of K\ and K 2 onto the corresponding invariant 
subspaces, the transfer matrix is block diagonalized. The dimensions of the blocks are given 
in Table |. We note that using these projectors leaves a few degeneracies inside some blocks. 
We also note that some eigenvalues are independent of the value of K 2 . In Appendix [A] 
all these 2 N / 2+l ^-independent eigenvalues occuring only when the size N is a multiple of 
four, are determined analytically. 

Having discarded the degenerate states and the ^-independent states in each repre- 
sentation one can perform the RMT analysis. Fig. [I] shows the level spacing distribution 
for wi = w 2 = 1.4 of the representation labelled R = 4, C = of iV = 16 (see Table |). 
One finds roughly an exponential distribution which is expected for an integrable model. 
Our explanation for the deviation from the exponential is the following: for an integrable 
system we do not sort the states according to all their "quantum numbers" (doing that will 
leave us with blocks of size 1). We therefore treat together eigenvalues belonging to states 
having different symmetries, but these states are only approximately independent. We also 
have performed the same numerical analysis on all eigenvalues computed with formula (AT) 



and the result is even worse. This result is characteristic of the free-fermion nature of the 
problem (note that the parameter {3 has been found always negative on the free-fermion 
variety of the eight- vertex model 0). The same form of the spacing distribution is found 
for any value of the Boltzmann weights, even for negative Boltzmann weights. However, if 
one of the Boltzmann weight is 'too large' or 'too small' some entries of the transfer matrix 
become huge and this leads to numerical instabilities in the diagonalization. On the other 
hand, if the Boltzmann weights are 'too close' to unity (decoupling limit) many eigenvalues 
are almost-degenerate which leads to difficulties in the unfolding procedure. We note finally 
that the critical point does not manifest itself in any manner on the spacing distribution nor 
on the spectral rigidity A 3 . 

B. Two-Dimensional Ising Model in a Magnetic Field 

We now investigate the case where a magnetic field is turned on. The Hamiltonian 
becomes: 

N-l M-l 

PHnM k u k 2 ) = pn° NM (K u k 2 )-kY, E <n* > ( 13 ) 

K being the field times the inverse temperature, and the transfer matrix reads 

T^K^K^K) = [V(K 2 )B(K)] 1 ^ 2 H(K 1 )[V(K 2 )B(K)] 1 ^ 2 , (14) 

where H and V are the same matrices as defined by Eq. (|T0| ) and @ and B(K) is a diagonal 
matrix with entries 
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JV— 1 

[B{K)} a>a =U W 

i=0 



(15) 



where = ±1 according to the i bit of the binary representation of a, and w = e . The 
spin reversal symmetry does not hold any more. 

It is known that this model is not Yang-Baxter integrable and its partition function has 
not yet been summed up. We present in Fig. ^| the spacing distribution and the spectral 
rigidity for a typical large representation (labelled R = 4 in Table |[) for L = 14. The 
temperature is T = 1.4 and the magnetic field is H = 0.8. The statistics is taken over 
1100 spacings. The spacing distribution clearly coincides with the Wigner Surmise. The 
agreement of the spectral rigidity A 3 and the number variance S 2 with same quantities 
computed for the GOE matrices is surprisingly good: it holds up to a value of L = 25 for 
the spectral rigidity and only up to L = 7 for the number variance, a value much larger 
than for other models of statistical mechanics 0. To appreciate how the magnetic field 
influences the spacing distribution we have plotted in Fig. | the best fitted value (3 of the 
Brody distribution Pp Eq. @ as a function of the Boltzmann weight associated to the field. 
The behavior of (3 is unambiguous: the magnetic field induces a Wigner type distribution 
of the spacing distribution. The drop of (3 is sharp as the magnetic field goes to zero. On 
Fig. m the negative value of the parameter (3 in zero field is due to the free-fermion nature 
of the two dimensional Ising model in absence of magnetic field. We interpret the small 
width of the peak as a size effect, and we claim that in the thermodynamic limit (3 is strictly 
zero only for a zero magnetic field. This underlines the very singular nature of integrability. 
Probing the value of (3 in the physical region of the parameter space, we did not have find 
other integrable points than the points where the magnetic field is zero. 



C. Three-Dimensional Ising Model 

We now investigate another archetypal case of non-integrability in statistical mechanics: 
the three-dimensional Ising model. This model can still be mapped onto a fermion problem 
like in two dimensions, but in three dimensions the fermions are correlated |20| and the 
partition function cannot be summed up in a closed form. 

We build the lattice by adding square isotropic layers. Be K 2 the interaction in the 
two directions inside the layers and K\ the interaction in the direction perpendicular to the 
layers (i.e. the cubic lattice is built in the direction of K\). We again have T^(Ki, K 2 ) = 
V 1 / 2 (K 2 )H(K 1 )V 1 / 2 (K 2 ). The matrix V(K 2 ) is a 2^ x 2 N * diagonal matrix with entries 

N-1 N-1 

[v(K 2 )] a>a = n n w ^ +i+a ^ } (i6) 

i=0 j=0 

where a^j = ±1 (the sums in the indices are taken modulo N). The matrix H(K{) is a 
symmetric matrix with entries 

N-1 

[H{Ki)U = II <" Ai • (17) 

i,j=0 
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We again need to project the transfer matrix into the parameter independent invariant 
subspaces. Without magnetic field we still have the spin reversal symmetry, but the space 
group is more involved than for the two-dimensional model, since the automorphy group of a 
square lattice is larger than the automorphy group of a ring (see Sec. ||). Using rectangular 
N x N' layers would have led to much simpler calculations since the symmetry group would 
have been simply Pjv ®Dn'- However, the size reduction of the matrix would have been less 
and consequently also the numerically accesible lattice sizes. It is worth noting at this point 
that the size effects are not controlled by the dimension of the subblocks of the transfer 
matrix, but by the size of the lattice. Moreover, we have studied the particular case iV = 4 
for which the generic symmetry group of the isotropic square lattice C^ v has to be replaced 



by a larger group (for the technical details see [lql). In contrast with the two-dimensional 
case, projecting onto invariant subspaces lifts all degeneracies within each subblock (note 
that using the generic C iv group would have left degeneracies inside some blocks). Table |lj 
shows the size of the invariant subspaces. 

The results of the RMT analysis are presented in Fig. f|. The eigenvalue spacing dis- 
tribution for the representations labelled R = 17 and R — 18 in Table are averaged 
together. After discarding some eigenvalues close to the edge of the spectrum we are left 
with 2000 spacings. We observe that the level spacing distribution is very close to a Wigner 
distribution. Also the spectral rigidity and the number variance are in agreement with the 
corresponding quantities of GOE matrices. We have paid special attention to the critical 
point, since it is an example of a nonintegrable, nevertheless critical, point. No special 
behavior is found. This strongly suggests that the statistical properties of the eigenvalues 
of transfer matrices are governed by the status of the model with respect to integrability 
rather than to criticality. 

To go from a two-dimensional model to a three-dimensional model we have recorded the 
value of (3 when the interaction w\ in the third dimension is turned on continuously. Fig. |5| 
summarizes the results: a vanishingly small coupling in the third direction induces level 
repulsion. It is remarkable that with such a small size N = 4 we already have a very abrupt 
variation of (3. This once again stresses the singular nature of integrability. 



D. Three-State Potts Model on a Square Lattice. 



We now turn to the case of the Potts model (see Ref. PI] and reference therein). This 
spin model is a generalization of the Ising model where the spins can take more than two 
values. The Hamiltonian is: 

N-l AI-l 

PH% M (K U K 2 ) = -Y,H {KA°i,v <n*+i) + K 2 5(a hJ , a i+1>j )) , (18) 

where 5 is a Kronecker symbol, (3 is the inverse temperature, K\ and K 2 are the coupling 
constants divided by the temperature and a G Z q can take q values. It has been investigated 
by many authors but its full solution is still a challenge. Using a duality relation one can 
localize a phase transition at the temperature T c = 1/ ln(l + Jq) for any number q of states 
|22|| . The model can be mapped onto a staggered six- vertex model, the parameters of which 



depend on Ki, K 2 and q pi. This six- vertex model is in general not integrable since there 
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are two different sets of Boltzmann weights, one for each sublattice of the square lattice. 
However, for a special combination of the parameters the two sets of Boltzmann weights 
are the same and, consequently, the partition function of the model can be calculated. This 
line where the partition function can be calculated turns out to be the critical line. Here, 
in contrast with the three-dimensional Ising model, the two notions of integrability and 
criticality coincide. 

We have numerically investigated the case of the three-state Potts model. For q = 3 
this model presents a second order phase transition which is not of the same universality 
class as the transition of the Ising model. For the isotropic case the transition is given by 
(e K — I) 2 = q = 3. It would have been interesting to study higher values of q for which the 
transition becomes first order, but the size of the transfer matrix is an exponential function 
of the number of states q, and values of q larger than 3 lead to extremely large matrices 
even for a small lattice size. The transfer matrix T has the same form (Eq. flTTD) as for the 
two-dimensional Ising model in absence of magnetic field. The matrix V(K 2 ) is a 3 N x 3 N 
diagonal matrix with entries 

N-l 

\y(K*)fo = n ^ K,Q!+i) , (19) 

i=0 

where aj = 0, 1 or 2 is the value of the i th spins in the spin configuration labelled by a. The 
matrix H(K\) is a symmetric matrix: 



H{K l ) = h{K 1 ) m , h(K 1 )= 1 Wl 1 I . (20) 



The space symmetry group is T*n as in the case of the Ising model, but the color symmetry 
group is S3 rather than S2 for the Ising model. The size reduction is better in the Potts 
case but far less than the exponential increase due to the three states. Tab. p II] gives the 
size of the blocks in the different invariant subspaces. On Fig. |6| we present two typical level 
spacing distributions. The upper one (a) is obtained at the critical value of the Boltzmann 
weight w* = 1 + v3 and the other one is obtained at a different value w = 1.4 far from the 
transition. It is obvious that at w* the distribution is very close to an exponential while 
for w w* this distribution is close to the Wigner surmise. Fig. [7] shows the rigidity A 3 
for the same values of the Boltzmann weights. We observe the same coincidence with the 
theoretical behavior expected for independent numbers and spectra of GOE matrices. The 
agreement with the GOE behavior extend up to a value L ~ 4. This value is much less than 
for the two-dimensional Ising model in a field, but are comparable with values obtained in 




some quantum models |14[ . We present on Fig. [g the behavior of the best fitted value (3 
of the parametrized distribution as a function of temperature. Several sizes N = 8, 9, 10 
and 11 are plotted. We first note that the spacing distribution has properly 'detected' the 
integrable point w*, which corresponds precisely to the minimum in the /3(w) curve. We 
also observe size effects: for larger size the downward peak is sharper than for smaller size. 
This suggests that in the thermodynamic limit the spacing distribution is a Poisson law only 
at the critical temperature. 
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We have again found that integrability leads to independent eigenvalues, whereas non- 
integrability leads to eigenvalue repulsion properly described by the spectral statistical prop- 
erties of the GOE: the analysis of the statistical properties of the spectrum of the transfer 
matrix can be used to find integrable points. 



E. Three-State Potts Model for Negative Boltzmann Weights 

We now look at the isotropic three-state Potts model with w < 0. The transfer matrix 
is still symmetric, but some entries become imaginary, since we have half-integer powers of 
negative numbers (see Eqs. (p!9|) and (|TT|)). We have numerically found that the spectrum 
is complex only when —2 < w < 0. So when w < —2 one can apply the standard RMT 
analysis. We have found a GOE spacing distribution as expected. When — 2 < w < the 
spectrum contains mostly complex conjugate eigenvalues. By spacing we now mean shortest 
euclidian distance between eigenvalues in the complex plane, and we use another unfolding 
procedure ||, Chap. 8.6]. We want to know if the unfolded eigenvalues are 'independent' or 
if they repel each other. The spacing distribution of independent points in d dimensions is 
easy to evaluate: for N points taken randomly on a d- dimensional hyper-sphere of radius 
N 1 ^, the probability that the distance between a point and its closest neighbor is s, is the 
probability that exactly one point is found at the distance between s and s + ds, and that 
the other N — 2 points are farther away than s: 



P/v(s)cfe oc s c 



1 - 



ds . (21) 



Taking the limit N — > oo one gets P(s) oc s d_1 exp (— Cs d ) {C is some constant). For d = 2 
one recovers precisely the Wigner law which is a well known distribution in mathematical 
statistics. In other words a Wigner law for a one- dimensional set of points means repulsion 
between these points, whereas for points in the plane the same Wigner law means indepen- 
dence. The probability that random points in the plane are close to each other is already 
small and a supplementary repulsion due to correlations of non-integrability will have less 
influence than in the case of a real spectrum. For nonsymmetric random matrices the joint 



probability distribution of the eigenvalues has been studied [23[] . Eigenvalue repulsion is still 
present, but a closed expression for the degree of repulsion is not known. 

On Fig. [|(b) we present a distribution of eigenvalue spacings for w = —1.5: the repulsion 
is clearly seen, since P(s) near the origin is smaller than for the Wigner law. We note that 
this distribution is not close to the eigenvalue spacing of the spectra of GUE matrices. We 
also have investigated another special point in the regime —2 < w < 0. The self-duality 
equation (wsd ~ I) 2 — Q h as a solution wsd = 1 — v^3 which lies in this regime. On 
Fig. |(a) we present the distribution of eigenvalues spacings at this value wsd = 1 — \/3. 
The agreement with the Wigner law, which here means independence, is quite good. This 
is expected since the point wsd is also integrable. 
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IV. DISCUSSION AND CONCLUSION 



We have numerically shown that the eigenvalues of the transfer matrix of an integrable 
spin model have many features of random independent numbers. On the contrary, the 
spectrum of the transfer matrix of a nonintegrable spin model has many features in common 
with the spectrum of a GOE matrix. In particular, the spacing distribution is an exponential 
law for integrable systems while it is close to the Wigner surmise for nonintegrable systems. 
Quantities involving more than two closest eigenvalues, like the rigidity, also show a GOE 
behavior on a quite large scale involving up to 25 eigenvalues in the case of the square 
lattice Ising model with field. The independence of eigenvalues for integrable models has 
been checked on the Ising model in two dimensions in absence of a magnetic field, as well as 
on the critical point of the three-state Potts model. We also have studied the transition point 
of three-dimensional Ising model, for which we have found eigenvalue repulsion and a good 
agreement of the spacing distribution with the Wigner surmise. This stresses the difference 
between criticality and integrability: the eigenvalue statistics is sensible to integrability and 
not to criticality. 

Using the eigenvalue statistics as a criterion of integrablity we support the hypothesis 
that the two-dimensional Ising model in a field and the three-dimensional Ising model are not 
integrable. The integrable models appear as very singular and isolated in parameter space. 
This is clearly seen on curves showing the parameter (3 as a function of a Boltzmann weight, 
where (3 is close to unity almost everywhere, except for the particular values where the model 
is integrable. With the sizes numerically tractable the variation of /3 is abrupt, and the size 
behavior suggests that in the thermodynamic limit the statistics changes discontinuously. 
To clearify this point, a more detailed study of the size effects is needed. It would also be 
interesting, but difficult, to study the Potts model with a large number of states to have a 
first order phase transition point which is integrable. 

To use the criterion of eigenvalue spacing statistics with new models (for example chi- 
ral models) one needs to study complex spectra in many cases. The distinction between 
independent eigenvalues for integrable models and repelling eigenvalues for nonintegrable 
models still holds. However, a repulsion between eigenvalues in one dimension is much 
easier to quantify than in two dimensions. Intuitively this can be understood since in two 
dimensions the eigenvalues are not restricted to a line and have naturally more space to 
avoid each other. The repulsion between eigenvalues in two dimensions has less effect than 
in one dimension. We have found that for the three-state Potts model and for a Boltzmann 
weight — 2 < w < the spectrum is complex and that the eigenvalue spacing distribution 
is characteristic of eigenvalue repulsion. However, when w is close to the negative self-dual 
value we recover a distribution close to the Wigner law indicating here eigenvalue indepen- 
dence associated to integrability. The numerical difficulties arising with complex spectra are 
of importance and we need to refine further the analysis, and especially the unfolding, to 
study convicingly chiral models. This is in progress. 

In these two papers we have shown numerically that the statistical properties of transfer 
matrix spectra of classical statistical mechanics models are related to the integrability of the 
model, extending hereby the field of application of RMT. This can be useful in the search 
for new integrable models. 
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APPENDIX A: 



It has been established in jn| that the spectrum of the row-to-row anisotropic Ising 



i N-l 



model for an even size N is given by: 

A r (Jfi, K 2 ) = [2 sinh (2K 1 )} N/2 exp I ~ £ r q e q J (Al) 

where r q = ±1 according to the q th digit of the binary representation of r (0 < r < 2^), 
and for all q with < q < N: 

cosh e q = cosh 2K 2 cosh 2K{ - sinh 2K 2 sinh 2K\ cos Q , ( A2) 

where Q = (2q + 1)tt/N if there is an even number of 1 in the base two representation of 
t and Q = 2qir/N otherwise. For e q , the positive root of Eq. (|A2|) has to be taken except 
for Q = where eo = 2(K* — K 2 ) and for Q = it where e^/ 2 = 2{K\ + K 2 ). K{ is the dual 
coupling constant defined by tanhA'i = e~ 2Kl . The derivative of A T (Ki, K 2 ) with respect 
to K 2 is 

dAAK u K 2 ) I ^ de q 

dK 2 = ~2 AAKl ' K2) g ^ (A3) 

with 

de q ^ sm h 2-^2 cosh 2K* — cosh 2A" 2 sinh 2K* cos Q , . . 

dK 2 sinhe 9 

It is easy to check that for an even number of particles: 

e q = ejv -i_ 9 and = ^-=i (0 < q < N/2) , (A5) 

while for an odd number of particles: 

e q = e N . q and ^ = ^ (1 < q < N/2) , (A6) 



and 



<9f 

eo = 2(*J - JT 2 ) , ^ = 2, (A7) 
e N/2 = 2{K{ + K 2 ) , -ftp- = -2. (A8) 
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When the size N is a multiple of 4 and when there are N/2 particles, then the number 
of particles is even and we must use the formula (|A5|). There are 2 7V//2 choices of r such that 
r q = — Tjv-i-g for < q < N/2. In that case one has: 



JV-l q N/2-1 g 

E^= E (r q -r N ^ q )^ = 0. (A9) 

g=0 0K 2 q=0 OK 2 



This gives 2 N ' 2 -^-independent states. 

N being still a multiple of 4 and taking N/2± 1 particles we now must use formulae (|A6|) 
and ( |A7| ). There are 2 x 2 N I 2 ~ X choices of r such that r q = —T^^ q for 1 < q < N/2 and 
tq = In that case one has 



N-l a N/2-1 o 

E^ = ±(^o + ^ /2 ) + £ (^-^-i-«)^ = o. (A10) 

q=0 0K 2 g=l ^2 



Collecting the three cases mentioned above we get 2 N ^ 2+1 ^-independent states. It is 
then simple to verify that the corresponding eigenvalues are 

2 N {smhK 1 ) p {coshK 2 ) N - p , (All) 

where p = N/2 — 1, N/2 or N/2 + 1 is the number of particles. 
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TABLE I. Two-dimensional Ising model: the symmetries of a periodic chain combined with the 
Z2 color symmetry. The dimensions an and degeneracies Ir of the invariant subspaces for N = 14 
and N = 16. R is an arbitrary label of the representations of the dihedral group, exp(i/c) and 
A are the eigenvalues of the corresponding translation and reflection operators (* means that the 
corresponding representation is not stable under the action of the reflection operator). The column 
an : c=o (resp. cir : c=i) refers to states which are even (resp. odd) under spin reversal. The column 
cir is the sum of the two preceding for the case where a magnetic field breaks the spin reversal 
symmetry. 
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Ising Model on a Periodic Cubic Lattice 








N=4 


dim=65536 




R 


Ir 




a>R,c=i 







1 


222 


180 


402 


1 


1 


50 


44 


94 


2 


1 


169 


180 


349 


3 


1 


33 


44 


77 


4 


2 


191 


192 


383 


•5 


2 


211 


192 


403 


6 


3 


186 


180 


366 


7 


3 


291 


300 


591 


8 


3 


183 


180 


363 


9 


3 


354 


300 


654 


10 


4 


460 


480 


940 


11 


4 


236 


224 


460 


12 


4 


460 


480 


940 


13 


4 


236 


224 


460 


14 


6 


397 


416 


813 


15 


6 


507 


480 


987 


16 


6 


447 


480 


927 


17 


6 


681 


672 


1353 


18 


8 


668 


672 


1340 


19 


8 


668 


672 


1340 


TABLE II. 


Three-dimensional Ising model: 


the symmetry of an isotropic square lattice 


with 


periodic boundary conditions combined with the Zi color symmetry. The notations are the 


same 


as in Table |. 














Three-State Potts Model 


on a Periodic Square Lattice 





N=ll dim=177147 

R k A Ir aR,C=0 Qfi,C=l «R,C=2 or 

1 1 1464 1342 2806 8418 

1 vr -1 1 1221 1342 2563 7689 

2 vr/11 * 2 2684 2684 5368 16104 

3 2vr/ll * 2 2684 2684 5368 16104 

4 37r/ll * 2 2684 2684 5368 16104 

5 4tt/11 * 2 2684 2684 5368 16104 

6 5tt/11 * 2 2684 2684 5368 16104 



TABLE III. Three-state Potts model: the symmetries of a periodic chain combined with the S3 
color symmetry. Same notations as in Table | but C has three possible values C = 0, 1 or 2. There 
is an extra twofold degeneracy for C = 2. 
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Isotropic 2D Ising Model without Field 

"I t 

1 - 
0.8 

2 0.6- h \ 



0.5 1 1.5 2 2.5 3 3.5 4 

s 

FIG. 1. The eigenvalue spacing distribu- 
tion for the row-to-row transfer matrix of the 
isotropic two-dimensional Ising model. The lin- 
ear size is N = 16 and the Boltzmann weight is 
w± = W2 = 1.4. The spacings are taken in the 
largest representation. The average is over 1400 
observations. 



Isotropic 2D Ising Model with Field 




FIG. 2. The eigenvalue spacing distribution 
P(s), the rigidity A3 and the number variance 
T? for the row-to-row transfer matrix of the 
isotropic two-dimensional Ising model in a field. 
The linear size is N = 14, the temperature is 
T = 1.4 and the magnetic field is H = 0.8. The 
spacings are taken in the representation labelled 
R = 4. The average is over 1100 observations. 
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Isotropic 2D Ising Model with Field 




0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 
w=exp(H/T) 

FIG. 3. The parameter f3 for the isotropic 
square lattice Ising model as a function of the 
Boltzmann weight w = exp(H/T) associated 
with the field. The other Boltzmann weights are 
kept constant w\ = W2 = 1.6 which determines 
the temperature as T ~ 2.1. Average over the 
representations R = and R = 4 for iV = 14 
(approx 1700 spacings). 

Isotropic 3D Ising model without Field 

1 ; 
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FIG. 4. The eigenvalue spacing distribu- 
tion for the row-to-row transfer matrix of the 
isotropic three-dimensional Ising model at the 
critical point. The linear size is N = 4. The 
spacings are averaged over the representations 
labelled 17 and 18. The average is over 2000 
observations. 



Anisotropic 3D Ising Model without Field 
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expfK,) 

FIG. 5. The best fitted [3 parameter as a fucn- 
tion of e Kl for the anisotropic three-dimensional 
Ising model. The average runs over 2000 spac- 
ings. The point e Kl = 1 corresponds to a 
two-dimensional model. 
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FIG. 6. Eigenvalue spacing distribution for 
the three-state Potts model (a) precisely at the 
critical point w\ = u>2 = 1 + v3 and (b) far 
from the critical point at w\ = W2 = 1.4. The 
data are obtained for N = 11 and the number of 
spacings is 2500 (a) and 1400 (b). 



Potts q=3 



Potts q=3 




Potts q=3 



J 
10 



■ (b) 










/ f » + 















2 4 6 
N 



10 



FIG. 7. The rigidity A3(L) for the same pa- 
rameters as in figure [f| 

Potts q=3 
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FIG. 8. The parameter (3 as a function of the 
temperature for a three-state Potts model for 
different lattice size. The minimum is at the 
critical point which is integrable. 
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FIG. 9. Two closest-distance distributions for 
the isotropic three-state Potts model with neg- 
ative Boltzmann weight: (a) for the self-dual 
Boltzman weight w\ = W2 = wsd = 1 — V?> 
and (b) for a different value w\ = W2 = —1.5. 
The average runs over about 2000 spacings. 
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